Electronic, vibrational and magnetic properties of a novel C48N12 aza-fullerene 
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We study the structural, electronic, vibrational and magnetic properties of a novel C48N12 aza- 
fullerene using density functional theory and restricted Hartree-Fock theory. Optimized geometries 
and total energy of this fuUerene have been calculated. We find that for C48N12 the total ground state 
energy is about -67617 eV, the HOMO-LUMO gap is about 1.9 eV, five strong IR spectral lines are 
located at the vibrational frequencies, 461.5 cm~^, 568.4 cm~^, 579.3 cm~^, 1236.1 cm~^, 1338.9 
cm~^, the Raman scattering activities and depolarization ratios are zero, and 10 NMR spectral 
signals are predicted. Calculations of diamagnetic shielding factor, static dipole polarizabilities 
and hyperpolarizabilities of C48N12 are performed and discussed. Our results suggest that C48N12 
may have potential applications as semiconductor components and possible building materials for 
nanometer electronics, photonic devices and diamagnetic superconductors. 



Since the discovery of Cgo [1], doped fuUerenes have 
attracted a great deal of interest in the research commu- 
nity of physics, chemistry and material engineering due to 
their remarkable structural, electronic, optical and mag- 
netic properties [2-5]. For example, it has been shown 
that the doped fuUerenes can exhibit large third-order op- 
tical nonlinearities and be ideal candidates as photonic 
devices including all-optical switching, data processing 
and eye and sensor protection [4]. Another example 
is alkali-doped Ceo crystals, which can become super- 
conductors [6,7], for example, at a critical temperature 
Tc = 30 K [6]. Besides the alkali metal doping, there 
is also another type of doping, named substitute dop- 
ing [2-5], i.e., substituting one or more carbon atoms of 
fuUerenes by other atoms. Over the past 10 years, boron 
and nitrogen atoms have been successfully used to re- 
place carbon atoms of Cgo and synthesize many kinds 
of doped Ceo molecules, Ceo-m-nXmYn [2-5,8]. A very 
efficient method of synthesizing C59N has been recently 
developed by Hummelen et al. [9]. Their method has 
led to a number of detailed studies of the physical and 
chemical properties of C59N [3-5], which suggests that 
the new type of substituted fuUerenes may be used as 
semiconductor components and possible building mate- 
rials for future nanometer electronics since their band 
gaps and electronic polarizations can vary with different 
substitute doping [2-5]. 

Very recently. Ceo with more than one nitrogen atom 
replacing a carbon atom in the cage has been synthesized 
by Hultman et al. [10]. They reported the existence of a 
novel C48N12 aza-fullerene [10,11], which has one nitro- 
gen atom per pentagon and the symmetry of the Sq point 
group. Hence, it would be interesting and significant to 
investigate and predict the structural, electronic, vibra- 
tional and magnetic properties of this new aza-fullerene 
from the viewpoint of practical applications. This forms 
the purpose of the present paper. In the following, we 



report ab initio calculations of the optimized geometric 
structure, electronic properties, vibrational frequencies, 
IR spectra, NMR shielding tensors, diamagnetic shield- 
ing factor, static dipole polarizabilities and hyperpolar- 
izabilities of the stable C48N12 aza-fullerene. Our results 
suggest that C48N12 may have potential applications as 
semiconductor components and possible building materi- 
als for nanometer electronics and diamagnetic supercon- 
ductors. 
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FIG.l: C48N12 geometry structure obtained 
from ab initio B3LYP geometry optimizations 
with a 6-31G basis set. The center numbers {5, 
9, 14, 21, 26, 30, 35, 39, 45, 50, 55, 60} are for 
nitrogen atoms and the others for carbon atoms. 
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Full geometry optimization and total energy calcula- 
tion of C48N12 were performed by using the NWChcm 
program [12,13], where we have employed the B3LYP [14] 
hybrid density functional theory (DFT) method that in- 
cludes a mixture of Hartree-Fock (exact) exchange, Slater 
local exchange [15], Becke 88 non-local exchange [16], 
the VWN III local exchange-correlation functional [17] 
and the LYP correlation functional [18]. Throughout the 
split- valence 6-31 G basis set is used in our calculations. 
All results are calculated to machine precision. As will 
be discussed later, the best total energy calculations with 
other basis sets differ from the 6-31G basis set results by 
less than one percent. This should define the accuracy of 



our calculations. 

The optimized geometric structure of C48N12 is pre- 
sented in Fig.l. As shown by Hultman et al. [10], the ab 
initio calculations show that there is only one nitrogen 
atom per pentagon and two nitrogen atoms preferentially 
sit together in one hexagon. The symmetry of the stable 
isomer of C48N12 is the Sq point group [11]. The radial 
distance from the zth atom to the center (i.e., the ori- 
gin) of C48N12 is listed in Table I. We find that there 
are 10 different radii, which suggest that C48N12 is an 
ellipsoid structure, different from Cgo where each carbon 
atom has equal radial distance from the center. 



Table I: Calculated net MuUiken charges Qi ( g = 1.602 x 10~^^), electrostatic potential Vi ( t; = 
2.7209 X 10^), electron density pi {rj = 1.08 x 10^^)) and radial distance Ri at the atom center number Uj 
in C48N12 by using B3LYP hybrid DFT method with a split-valence 6-31G basis set. 



Hi 


Atom 


Ri 


Qi/q 




Pi/V 






[ nm] 


[C] 


[V] 


[ C/m^ ] 


{ 8, 20, 22, 44, 46, 56} 


C 


0.3495 


0.262 


-14.6663 


117.999 


{ 10, 17, 25, 41, 49, 59} 


C 


0.3508 


0.255 


-14.6674 


118.029 


{ 4, 15, 27, 34, 40, 54} 


C 


0.3541 


0.249 


-14.6679 


118.027 


{1, 13, 16, 31, 38, 51} 


c 


0.3535 


0.215 


-14.6686 


118.058 


{ 7, 19, 23, 43, 47, 57} 


c 


0.3448 


0.210 


-14.6631 


118.004 


{ 6, 18, 24, 42, 48, 58} 


c 


0.3470 


0.199 


-14.6694 


117.992 


{ 2, 12, 29, 32, 37, 52} 


c 


0.3545 


0.002 


-14.7041 


118.109 


{ 3, 11, 28, 33, 36, 53} 


c 


0.3549 


-0.017 


-14.7026 


118.103 


{ 5, 14, 30, 35, 39, 55} 


N 


0.3595 


-0.669 


-18.2639 


191.331 


{9, 21, 26, 45, 50, 60} 


N 


0.3531 


-0.705 


-18.2547 


191.325 


Carbon atoms in Cgo 


c 


0.3535 





-14.7146 


118.181 



The net MuUiken charges of carbon and nitrogen atoms 
in C48N12 are listed in Table I. It is found that there 
are two types of nitrogen atoms in the structure. The 
net MuUiken charges on the two types of nitrogen atoms 
arc —0.67 q and —0.71 q. The carbon atoms separate 
into two groups, one-fourth of the carbon atoms having 
net MuUiken charges in the range of —0.02 q to 0.01 q 
and the remaining three- fourth in the range of 0.19 q to 
0.27 q. Although the MuUiken analysis cannot estimate 
the atomic charges quantitatively, their signs can be es- 
timated [19]. From these results, we find that the doped 
nitrogen atom exists as an electron acceptor whose net 
MuUiken charge is negative, and 7/8 of the carbon atoms 
as electron donors whose net MuUiken charges are posi- 
tive. We should mention that we also performed calcula- 
tions of net MuUiken charges of carbon and boron atoms 
in C48B12. It is found that the doped boron atom exists 
as an electron donor and carbon atom as an electron ac- 
ceptor (more detailed results for C48B12 will be reported 
in another paper [20]), which is consistent with the exper- 
imental result [8]. Therefore, C48N12 and C48B12 have 
opposite electronic polarizations, while Cgo is isotropic. 
In the case of doping into silicon, we know that the doped 
phosphorous (the V family in the periodic table) exists 
as a donor, while the doped boron (the III family) exists 



as an electron acceptor. Thus, the results for C48N12 and 
C48B12 differ greatly from that for silicon. We propose 
that C48N12 and C48B12 may be important components 
of semiconductors with opposite electronic polarizations. 

The optimized C-C and C-N bond lengths in C48N12 
are listed in Table II. We find that there are six different 
C-N bond lengths and nine different C-C bond lengths, 
which are in agreement with the calculations of Stafstrom 
et al. [11]. The calculated C-C bond length s in C50 are 
also shown in Table II. It is seen that Ceo has only two 
kinds of bond lengths, i.e., one single C-C bond with 
0.14445 nm and one double bond with 0.13944 nm, which 
are in good agreement with the results measured by neu- 
tron scattering (0.145 nm and 0.1391 nm) [21] or nuclear 
magnetic resonance ( 0.146 nm and 0.140 nm ) [22]. It 
is obvious that most of C-C bond lengths in C48N12 , due 
to the electronic polarizations, are different from those in 
Ceo- 

The distribution of bond angles (C-C-C, C-N-C, C-C- 
N) in C48N12 is shown in Fig.2. The C-C-C and C-N-C 
(and C-C-N) bond angles separate into three groups: one 
in the range of 107° to 109°, one in the range of 116° to 
119° and the other one in the range of 119.5° to 122°. In 
contrast, there are only two kinds of bond angles, 108° 
and 120°, in Ceo- Hence, only few bond angles in C48N12 
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are similar to those of Ceo- 

Table II: Calculated bond lengths in C48N12 
by using B3LYP hybrid DFT method with a 
split-valence 6-31G basis set. 



Bond Center Number 



Bond Length 



c- 


c 


(3, 4) (11, 15) (27, 28) 
(33, 34) (36, 40) (53, 54) 


0.13971 


c- 


c 


(6 , 10) (17, 18) (24, 25) 
(58, 59) (41, 42) (48, 49) 


0.14004 


c- 


c 


(7,8 ) (19, 20) (22, 23) 
(43, 44) (46, 47) (56, 57) 


0.14068 


C- 


c 


(1 , 2) (12, 13) (16, 29) 
(31, 32) (37, 38) (51, 52) 


0.14125 


C- 


c 


(1 , 52) (2 , 31) (12, 38) 
(13, 29) (16, 37) (32, 51) 


0.14216 


c- 


c 


(6 , 7) (18, 19) (23, 24) 
(42, 43) (47, 48) (57, 58) 


0.14217 


c- 


c 


( 3, 10) (11, 59) (17, 33) 

(25, 36) (28, 49) (41, 53) 


0.14346 


C- 


c 


( 4, 46) ( 8, 15) (20, 40) 
(22, 54) (27, 44) (34, 56) 


0.14354 


c- 


c 


( 2, 3) (11, 12) (28, 29) 
(32, 33) (36, 37) (52, 53) 


0.14488 


N- 


c 


(7, 60) (9, 47) (19, 26) 
(21, 57) (23, 45) (43, 50) 


0.14189 


N- 


c 


(8 , 9) (20, 21) (22, 26) 
(44, 45) (46, 50) (56, 60) 


0.14149 


N- 


c 


(9 , 10) (17, 21) (25, 26) 
(41, 45) (49, 50) (59, 60) 


0.14226 


N- 


c 


(4 , 5) (14, 15) (27, 30) 
(39, 40) (34, 35) (54, 55) 


0.14286 


N- 


c 


(1 , 5) (13, 14) (16, 30) 
(31, 35) (38, 39) (51, 55) 


0.14315 


N- 


c 


( 5, 42) ( 6, 35) (14, 48) 
(18, 55) (24, 30) (39, 58) 


0.14317 


C- 


c 


carbon in Ceo 


0.13944 


c- 


c 


carbon in Ceo 


0.14445 
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FIG.2: ab initio B3LYP/6-31G DFT calcu- 
lation of C-C-C (open circles) and C-N-C (or 



C-C-N) (filled circles) bond angles in C48N12. 
The center numbers are labeled in Fig.l. 

Table III: Total energy calculations of 
C48N12 aza-fuUerene and Ceo by using 
B3LYP hybrid DFT method with a split- 
valence 6-3 IG basis set. 





Ceo 


C48N12 


total DFT energy (eV) 


-62192.890 


-67617.316 


one electron energy (eV) 


-537319.738 


-577022.480 


Coulomb energy (eV) 


256457.119 


274997.383 


exchange-correlation energy (eV) 


-8887.049 


-9361.705 


nuclear repulsion energy (eV) 


227556.778 


243769.848 


LUMO (eV) 


-3.382 


-2.707 


HOMO (eV) 


-6.216 


-4.633 


binding energy (eV/atom) 


8.676 


8.188 


first ionization energy (eV) 


7.548 


5.993 



The results of total energy calculations of C48N12 aza- 
fuUerene are summarized in Table III. For comparison, 
the results of Ceo are also listed in Table III. Their DFT 
orbital energies with the orbital symmetries are shown 
in Fig. 3. Because of the valency of the doped nitrogen 
atoms, the electronic properties of C48N12 are greatly 
different from Ceo- Since the icosahedral symmetry of 
Ceo is lost by the substituted doping, each energy level 
in C48N12 splits. In detail, we find that the ground state 
energy of C48N12 is -67617.316 eV, which is 5423 eV lower 
than that of Ceo. For Ceo, the gap obtained from the on- 
set of the occupied fivefold-degenerate /i„ energy level 
(i.e., HOMO) and the threefold-degenerate tiu level (i.e., 
LUMO) is 2.834 eV, which is in agreement with the ex- 
periment [26]. For C48N12, the HOMO is a doubly de- 
generate level of cig symmetry and the LUMO is a sin- 
gle level of Qu symmetry, and the corresponding LUMO- 
HOMO gap is 1.926 eV, which is 0.908 eV smaller than 
that of Ceo- These results indicate the possibility that 
C48N12 can be the components of semiconductors with 
small band gaps. Since C48N12 is isoelectronic with Cqq"^, 
it corresponds to a complete filling of the tiu and tig lev- 
els of Ceo [10]. From Fig.3, it is seen that the HOMO 
of C48N12 appears at about 1.583 eV above that of Ceo- 
This difference in the HOMO levels is in agreement with 
the calculation of Stafstrom et al. [11] and corresponds 
to the difference between the first ionization energies of 
C48N12 and Ceo- These results suggested that the aza- 
fulerene C48N12 is a very good electron donor [11]. The 
binding energies (i.e., the difference in the total energies 
of the whole system and the fragments) are also listed in 
Table III. We find that the binding energy for C48N12 is 
8.188 eV per atom, which is about 0.5 eV smaller than 
that (=8.676 eV per atom) of Ceo- Hence, as synthe- 
sized by Hultman et al. [10], C48N12 would be a stable 
aza-fuUerene. 
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Fig.3: B3LYP/6-31G DFT orbital energies 
Efc at the kth cigcnstate. Open circles with 
labeled orbital symmetries are for Ceo- Filled 
circles with orbital symmetry Ug and filled 
triangles with orbital symmetry a„ are for 
C48N12. LUMO and HOMO are shown. 



The interaction energy of a nucleus of magnetic mo- 
ment n with a weak external magnetic field B is given by 
E = —/X • B(l — a), where a is the diamagnctic shield- 
ing factor (DSF) [24,25] arising from induced magnetic 
fields at the nucleus. In Fig. 4, wo present the calculated 
DSF a of carbon and nitrogen atoms in C48N12 by using 
the NWchem program [12,13] with the B3LYP hybrid 
DFT method and the split- valence 6-3 IG basis set. In 
comparison with the calculated results of Cgo shown in 
Fig. 4, we find that the average DSF of carbon atoms 
in C48N12 is enhanced by about 3 %. This enhancement 
can be explained as follows. From Table I, we find that 
the electron density at each carbon atom in C48N12, in 
comparison with the Cgo's which is evenly distributed, 
is decreased due to the substitute doping. Correspond- 
ingly, as shown in Table I, the electrostatic potential at 
each carbon atom in C48N12 is also increased. The Lamb 
formula [24] shows that the DSF a is proportional to the 
electrostatic potential V produced at the nucleus by the 
electrons. Hence, the DSF a of carbon atoms in C48N12 is 
enhanced. Our calculated results suggest that if C48N12 
is doped with alkali metal [6,7] and used as a building 
compound for a superconductor material, then it may be 
a good diamagnetic superconductor because of the en- 
hanced DSF in the carbon atom. 



Fig. 4: Calculated diamagnetic shielding fac- 
tor a of carbon (open circles) and nitrogen 

(filled circles) atoms in both Cgo and C48N12. 
The center numbers are labeled in Fig.l. 



Next, we performed calculations of vibrational frequen- 
cies, IR spectrum , Raman scattering activities and de- 
polarization ratios, and NMR shielding tensors of C48N12 
by using the Gaussian 98 program [13,23] with B3LYP/6- 
31G hybrid DFT method. For C48N12, there are 3n- 
6 = 174 vibrational modes for n=60 (the number of 
atoms). The corresponding vibrational frequencies at 
298.15 K and 1 atmosphere of pressure arc shown in 
Fig. 5, which also presents the IR intensity at each vibra- 
tional frequency. Since experimental IR spectroscopic 
data do not directly indicate the specific type of nu- 
clear motion producing each IR peak, we do not give 
here the normal mode information for each vibrational 
frequency and the displacements of the nuclei corre- 
sponding to the normal mode. From Fig. 5, we find 
that C48N12 has five strong IR lines with the intensities 
8.44 X 10" m/mol, 6.81 x lO" m/mol, 5.93 x 10* m/mol, 
5.15 X 10* m/mol and 4.2 x 10* m/mol, being located 
at the vibrational frequencies 461.5 cm^^, 1236.1 cm^^, 
1338.9 cm~^, 579.3 cm~^ and 568.4 cm~^, respectively. 
The vibrational frequencies, 461.5 cm~^, 1236.1 cm~^, 
1338.9 cm"-'^, and 568.4 cm"-*^, correspond to two degen- 
erate normal modes. Also, it is found that the zero-point 
correction to the total electronic energy is about 10.1 
eV, and the thermal corrections to the total electronic 
energy, enthalpy, Gibbs free energy are 10.7 eV, 10.8 eV 
and 9.0 eV, respectively. By contrast, the calculated re- 
sult for Ceo is shown in the inset in Fig. 5, where we ob- 
serve four strong IR spectral lines with the intensities 
2.76 X 10* m/mol, 8.2 x 10^ m/mol, 1.0 x 10* m/mol and 
1.72 X 10* m/mol at the vibrational frequencies 553.0 
cm-\ 586.5 cm-\ 1208.5 cm'^ and 1479.3 cm-\ re- 
spectively, which are in good agreement with the experi- 
ment (528 cm-\ 577 cm'^, 1183 cm"! and 1429 cm"!) 
measured by Kratschmer et al. [27]. 
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FIG. 5: ab initio B3LYP DFT calculation 
of vibrational frequency and IR intensity of 
C48N12 with a split-valence 6-31G basis set. 
Inset is the result for Ceo calculated with 
B3LYP/6-31G. 

Table IV: Static dipole polarizabilities a ( 
1 au = 0.5292 x IQ-^o ni3) of C48N12 and Ceo 
by using B3LYP hybrid DFT method with a 
split-valence 6-31G basis set. 



FuUerene 








axz 


Otyz 


Q-zz 




[ au] 


[ au] 


[ au] 


[au] 


[au] 


[au] 


Ceo 


460.84 


0.00 


460.84 


0.00 


0.00 


460.84 


C48N12 


442.55 


0.00 


442.54 


0.00 


0.00 


445.88 



The static dipole polarizability measures the ability 

of the valence electrons to find an equilibrium configu- 
ration which screens a static external field [28]. Hence, 
molecules with many delocalized valence electrons should 
display large values of static dipole polarizabilities. In 
Table IV, we list the calculated static dipole polarizabil- 
ities of C48N12. For comparison, the results of Cgo are 
also listed in Table IV. Obviously, large values of static 
dipole polarizabilities are found for both C48N12 and Ceo- 
Actually, this is expected since the sp structure for car- 
bon chains is rich of tt bonds delocalized along the en- 
tire body of the system [28] . Our calculated polarizabil- 
ities for Cgo arc in agreement with those of Fowler et al. 
[29]. From Table IV, we notice that only axx, cxyy and 
Uzz are different from zero, i.e., there exist three mutu- 
ally perpendicular directions in the molecule for which 
the induced dipole moments are parallel to the electric 
field. The locus of points formed by plotting Xj^fa in 
any direction from the origin yields a surface called the 
polarizability ellipsoid [30] whose axes are x, y, and z. 
Obviously, because the electron distribution for Cgo, as 
discussed before, is isotropic in the x, y and z directions, 
there is no difference between z and x (or y) components 
of the polarizabilities a and the polarizability ellipsoid 
becomes a sphere. However, this difference for C48N12 is 
pronounced because of the electron polarization caused 



by the electron acceptor, N atom. The polarizability of 
C48N12 is the same in both x and y directions. Thus, the 
polarizability ellipsoid in C48N12 becomes a rotational 
ellipsoid with two equal axes x and y. 

It is known that the molecular vibration must be 
accompanied by a change in the polarizability of the 
molecule in order for a molecular vibration to be Ra- 
man active [30]. Although the polarizability ellipsoid 
may have higher symmetry than the molecule, all the 
symmetry elements possessed by the molecule will also 
be possessed by the ellipsoid [30]. Hence, if the polariz- 
ability ellipsoid is changed in size, shape, or orientation 
as a result of vibrational motion, a Raman spectrum will 
result [30]. In our case, the calculated Raman scattering 
activities and depolarization ratios for C48N12 are zero. 
Thus, its vibration is not Raman active, which implies 
that the polarizability of the molecule does not change 
in a molecular vibration, and a measurement of its depo- 
larization ratio cannot provide a means of distinguishing 
totally symmetrical vibrations from the rest. 

The first-order hyperpolarizability /? [31] of C48N12 is 
also calculated and found to be zero, the same as that of 
Ceo- This is actually expected since they display inver- 
sion symmetries. Consequently, this aza-fuUcrcne C48N12 
cannot produce the second-order nonlinear optical (NLO) 
interactions. Based on the available studies of the third- 
order optical nonhnearities of C59N and C58N2 [4], we 
expect that C48N12 may exhibit large second-order hy- 
perpolarizability 7 because of its large electronic polar- 
izations and small HOMO-LUMO gap. This suggests 
that C48N12 may be one of the molecules, for which we 
are looking to build third-order NLO materials for pho- 
tonic applications. 

Modern high-field, multipulse NMR spectroscopy has 
proven to be an exceptionally powerful technique in char- 
acterizing molecular systems and structures. Here we 
calculate the NMR shielding tensor of carbon and ni- 
trogen atoms in C48N12 aza-fullerene. Since no cur- 
rent functionals include a magnetic field dependence, the 
DFT methods do not provide systematically better NMR 
CS results than Hartree-Fock theory [23]. Hence, we 
compute the NMR shielding tensor of C48N12 with the 
gauge-including atomic orbital (GIAO) method [32] by 
using Gaussian 98 program [23] in the level of restricted 
Hartree-Fock (RHF) theory with the split-valence 6- 
31G basis set. Table V summarizes the results for the 
isotropic and anisotropic absolute shielding constants of 
the carbon and nitrogen atoms in C48N12. The isotropic 
shielding constant Ajgo is defined as Ai^o — {l^xx^ l^yy^ 
^zz)l'i [33], where Ay (i,j = x, y, z) is the component 
of the shielding tensor. The anisotropic shielding con- 
stant Aoniso, an indication of the quality of the shielding 
tensor, is defined as ^aniso = A3 — (Ai + A2)/2 [33], 
where Ai < A2 < A3 are the eigenvalues of the sym- 
merized shielding tensor. It is seen that the NMR shield- 
ing constants separate into eight groups for carbon atom 
and two groups for nitrogen atoms. Hence, we predict 
that there are 10 groups of spectral signals in the NMR 
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spectroscopy of C48N12. In contrast, the calculated ab- 
solute isotropic and anisotropic NMR shielding constants 
for Ceo are 54.10 ppni and 180.09 ppni, respectively, and 
there is only one NMR spectral signal. In order to com- 
pare the predicted values of Ceo to experimental results, 
we also list the calculated NMR shielding constants of 
carbon atoms in tetramethylsilane (TMS). We find that 
the isotropic NMR shielding constant of carbon atoms in 
Ceo, relative to that in TMS, is 146.7 ppm, which are 
in good agreement with the experimental result (=142.7 
ppm) of Ceo [34]. 

Table V: Restricted Hartree-Fock calcula- 
tions of the isotropic Ajgo and anisotropic 
^aniso absolute shielding constants of the 
carbon and nitrogen atoms with a split- 
valence 6-31G basis set in the GIAO method 
for C48N12 aza-fullerene, Ceo and tetram- 
ethylsilane (TMS). 



Center Number 


Atom 




^aniso 






ppm ] 


[ ppm ] 


{1, 13, 16, 31, 38, 51} 


C 


44.0 


173.0 


{4, 15, 27, 34, 40, 54} 


C 


57.9 


154.9 


{7, 19, 23, 43, 47, 57} 


C 


59.4 


152.2 


{10, 17, 25, 41, 49, 59} 


C 


59.5 


133.5 


{2, 12, 29, 32, 37, 52} 


c 


67.9 


169.2 


{8, 20, 22, 44, 46, 56} 


c 


70.7 


150.8 


{3, 11, 28, 33, 36, 53} 


c 


76.6 


141.2 


{6, 18, 24, 42, 48, 58} 


c 


78.2 


110.5 


{9, 21, 26, 45, 50, 60} 


N 


95.3 


202.6 


{ 5, 14, 30, 35, 39, 55} 


N 


114.4 


171.5 


all carbons in Ceo 


C 


54.1 


180.1 


all carbons in TMS 


C 


200.8 


21.4 



In this work, we only focus on the split-valence 6-31G 
basis set and discuss only the frequency analysis at 298.15 
K and 1 atmosphere of pressure. Detailed frequency cal- 
culations and thermochemical analysis at different tem- 
peratures and pressures, the basis set effects including 
minimal basis sets ST0-3G and polarized basis sets 6- 
31G* [23], and other theoretical methods, for example, 
MP2 [23] and CSGT method [33], wll be discussed in 
a full paper [35]. For example, the total electronic en- 
ergies of C48N12 with respect to basis sets ST0-3G, 3- 
21G, 6-31G and 6-31G* are -66795.835 eV, -67263.725 
eV, -67617.316 eV and -67637.659 eV, respectively, which 
shows that improving the basis set from ST0-3G to 6- 
31G* gives a significantly better total energy and the use 
of polarization functions is not cpsccially important for 
total energy calculations, i.e., 6-31G basis set can give a 
good estimate of the total energy of C48N12. 

In summary, we have performed ab initio calculation 
of optimized geometries, electronic properties, diamag- 
netic shielding factor, vibrational frequencies, IR spec- 
trum, NMR shielding tensor, dipolc polarizabilities and 
hyperpolarizabilities of the novel C48N12 aza-fullerene. 
Some calculated results of C48N12 are compared with 



those of Ceo at the level of the same theory. It is found 
that this novel aza-fullerene has some remarkable fea- 
tures which arc different from and can compete with 
Ceo- This aza-fullerene may have potential applications 
as semiconductor components and possible building ma- 
terials for nanometer electronics, photonic devices and 
diamagnetic superconductors since their band gaps are 
small, diamagnetic shielding factor in carbon atom can 
be enhanced and electronic polarizations vary largely due 
to the doping effect. 
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